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Abstract — In this paper we study the distribution of the scaled 
largest eigenvalue of complex Wishart matrices, which has diverse 
applications both in statistics and wireless communications. 
Exact expressions, valid for any matrix dimensions, have been 
derived for the probability density function and the cumulative 
distribution function. The derived results involve only finite sums 
of polynomials. These results are obtained by taking advantage 
of properties of the Mellin transform for products of independent 
random variables. 

Index Terms — Communication systems; performance analysis; 
eigenvalue statistics; the Mellin transform. 

I. Introduction 

Eigenvalue statistics of Wishart matrices play a key role in 
the performance analysis and design of various communication 
systems. Among these, the distribution of Scaled Largest 
Eigenvalue (SLE), defined as the ratio of the largest eigenvalue 
to the normalized sum of all eigenvalues, has been shown to 
be an important measure. The applicability of the SLE spans 
from classical problems in statistics |[T1, ||2l, IQ, Q to modern 
applications in wireless communications f5\, f6\, f7), fSl, f9\. 
Classical problems include testing the presence of interactions 
in a two-way model |jj and testing the equality of eigenvalues 
of certain matrices against various of alternatives (2], f3l, 
Q. Contemporary applications in wireless communications 
include non-parametric detection in array processing ||5] and 
spectrum sensing in cognitive radio networks |6|, Q, O, 
||9i . Specifically, for spectrum sensing applications, the SLE 
is formulated as a test statistics, which is first proposed 
by IS and further investigated in [71, E), M- The SLE based 
detector is the best known detector for single source detection, 
outperforming several classical detectors in reaUstic sensing 
scenarios |6l, fT^, f9|. Despite the importance of the knowl- 
edge of the SLE, existing results on its statistical properties are 
rather limited. In this paper, we aim to address this problem 
by deriving simple and exact expressions for the Probability 
Density Function (PDF) and Cumulative Distribution Function 
(CDF) of the SLE. 

The rest of this paper is organized as follows. In Section [U] 
we formally define the scaled largest eigenvalue of Wishart 
matrix followed by a concise survey on existing results. 



as well as the closed-form coefficients. Numerical examples 
are provided in |IV] to verify the derived results. Finally in 
Section [V] we conclude main results of this paper 

II. Definitions, Prior Results and Contributions 
Define a K x N (K < N) dimensional random matrix 
X with independent and identically distributed (i.i.d) complex 
Gaussian entries, each with zero mean and unit variance. The 
K X K Hermitian matri?([^ R = XX^^ follows a complex 
Wishart distribution with N degrees of freedom (d.o.f). We 
denote the ordered eigenvalues of R as Ai > A2 > ... > 
Xk > 0, and the normalized trace of R as T = tr{R}/iir — 
^ij 1-^- The scaled largest eigenvalue of R is formally 
defined as the ratio of its largest eigenvalue to its normalized 
trace, i.e., 

Ai Ai 

T' 



X := 



(1) 



Section III is devoted to deriving the exact SLE distribution 



where it can be verified that x € 

The distribution of X has been the subject of intense study 
in the literature. An exact expression for the distribution of X 
in terms of a high dimensional integral has been proposed 
in [l]. In |2|, a relation between Laplace transforms of 
random variables X and Ai was established. By symbolically 
inverting the Laplace transforms, some representations for 
the distribution of X were derived in |i3J, [i4J. Whilst being 
exact, these representations IT], (|2l, JSj, ^ can only be 
evaluated numerically for small values of K and N due to 
their unexplicit and complicated forms. Recently, motivated 
by its application in spectrum sensing, several asymptoticaj^ 
distributions of X have been derived lU, Q, ||9] via random 
matrix theory. Although these results are easy to compute, 
their accuracy can not be guaranteed for not-so-large K and 
N. As an example, in Fig. [T] we illustrate the accuracy of 
an asymptotic result based on Tracy- Widom distribution ('TW 
based approx.') from Q and an improved version ('TW based 
approx. with correction') from |8| with a typical choice of 
parameters in spectrum sensing K = 4 and N = 100. 

"The operator denotes conjugate-transpose. 

''Asymptotic in the sense that the matrix dimensions go to infinity while 
their ratio is kept fixed, i.e. K — > 00, N 00 and K/N — > r £ (0, 1). 

■^Corresponding to a situation of a sensing device with 4 antennas with 100 
samples per antenna. 
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Fig. 1. Accuracy of some existing asymptotic approximations to the SLE 
distribution: K = A, N = 100. Tlie 'TW based approx.' refers to the result 
of (7j and the 'TW based approx. with correction' refers to the result of [S]- 



From Fig. [T] we can see that the approximation errors of 
both asymptotic resuhs are non-negligible. Note that for other 
applications discussed in Section |lj the K, N values are often 
smaller than in spectrum sensing applications. For example, 
in |1| the authors considered K < A and < 100 in all the 
simulations. In addition, for the application considered in |5| 
it was remarked that choosing K > 2 does not result in any 
performance improvement and there N was always chosen to 
be no more than 20. Therefore, the approximation accuracy 
may become even lower when using the existing asymptotic 
results for the above applications. 

In this work, we derive exact expressions for the SLE 
distribution, valid for arbitrary K and N, as a finite sum 
of polynomials with some unknown coefficients. Closed-form 
expressions for the coefficients are obtained for the most useful 
system configurations of K = 2, 3, 4 with any N. These 
results are simple to calculate and do not involve any integral 
representations. To obtain these results, we adopt a novel 
approach based on the Mellin transform, which eliminates the 
need to handle correlations between Ai and T. The derived 
results yield a useful analytical tool in applications involving 
statistical properties of the SLE. 

III. The SLE Distribution 

A. Exact Distribution 

Although there exists intractable correlation between ran- 
dom variables Ai and T, it has been proved in f5\ that the 
random variables X and T are independent. As such, Ai equals 
the product of the independent random variables X and T. 
By this independence, the (z — l)th moment of Ai can be 
represented as 



E[xr 



E[{XTY 



E[X'~^]E[T' 



(2) 



Moreover, the {z — l)th moment of a random variable x, with 
PDF p{x), equals its Mellin transform as 



^p{x)dx -.^ Mz[v{x)], 



(3) 



where denotes the Mellin transform operation. Define 

fx^{x), frix) and fx{x) as the PDFs of Ai, T and X 
respectively. We have 



M.ihAx)] ^ MMx)]M.[fx(x)]. 



(4) 



By the Mellin inversion theorem, the PDF of X can be 
uniquely determined by the following contour integral 



fx{x) 



1 

27ri 



MMT{x)r'- 



(5) 



In principle, the above Mellin inversion integral can be eval- 
uated by using the residue theorem. Note that, following the 
above Mellin transform framework, a related distribution of the 
trace to the smallest eigenvalue has been derived recently 1 10|. 
The PDF of Ai admits the following representation [1 IJ , 

K {N+K)i-2i^ 



hAx) = J2' 



j = N-K 



(6) 



where Cij denotes the unknown coefficients. Closed-form 
coefficients formulas will be derived in the next subsection. 
Meanwhile, numerical algorithms are also available in ifTTl . 
(12] to calculate Cij for a given K and N. 

In order to apply the Mellin transform framework, we first 
need to calculate A4z[f\i{x)], which equals 



MAfM{x)] = Y. E 

1=1 j = N-K 



r(z + j) 



(7) 



It is well known that the sum of all eigenvalues of R, ^^j^ Aj, 
follows central Chi-square distribution with 2KN degrees of 
freedom, therefore the PDF of T = (^^i / K can be 
obtained as 



Mx) 



K 



KN 



{KN 



Its Mellin transform is 



MAhix)] = 



[KN - 1)! 



T{z + KN -I) 



(8) 



(9) 



Inserting (jTjl and (j9]) into the Mellin inversion integral (jSj) 
we have 



(KN ^\\ ^ (N+K)i-2^^ 



K 



i=l j=N-K 



-^A{x,z) (10) 



where 

A{x,z) 



1 
27ri 



T{z + KN-l) 



Az. (11) 
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' dz. 



(12) 



fx{x) 



{KN - 



{KN-j-2)\ 



(13) 



where 



By definition of the Meijer's G function ifTSll . as shown where 



in ( 12i on top of this page, the function A{x, z) can now be 



represented as 



A(a:,z) = G}:?( ^ 



KN - 1 
3 



By using the fact that 



-h-l 



(a- 6-1)! 

where 0{-) denotes the Heaviside step function 



(15) 



(16) 



X < 
X > 



(17) 



the PDF of X in ( 10 1 simplifies to the expression shown in ( 13 1 



on top of this page. 

We now focus on the CDF. By definition, the CDF of X, 
equals 

^-(2/) = SSE E *""--^c,,i3(y), (18) 



1=1 j=N-K 



where 



B{y) 



KN-j-2 



[l-f)dx 



{KN-3-2)\ 



(19) 



and y £ [l,oo). Using the definition of the hypergeometric 
function 2Fi{a,b; c; x) = J2^=o ^°('c)^!^" '^here (a)„ = 
r{a + n)/r{a) defines the Pochhammer symbol, the function 
B{y) becomes 

B{y) = C{y)e (^^-y^+C (^^) 9 (^y - - C{1), 



C{y) 



{KN~3-2)\{j + iy 
2F1 (-KN + 3+2,3 + 1- J + 2- 



V .(21) 



Since the parameters a, b and c of take integer values, the 
hypergeometric function can be simplified to a finite sum 
of polynomials. Consequently the function C(j/), after some 
manipulations, equals 

1 



Ciy) = I f 



KN-j~2 KN-j- 

E 

g=0 



i-i/KYij + q + l)- 
{KN - j - 2- q)\q\ 



-y 



(22) 

Inserting pO] ) into ( fTSj l, the CDF expression is summarized 
as ^T4\ on top of this page. 

Note that both the exact PDF expression ( pj) and CDF 



expression ( 14 1 involve the unknown coefficients Ci j inherited 



from (|6]l. Closed-form expressions for the coefficients Cij will 
be derived for K < 4 with arbitrary N. For other K values, 
one has to resort to numerical techniques ifTTl . lfT2l to obtain 
the values of Ci,j. 

B. Closed-form Coefficients 

In order to circumvent possible computational burden when 
using the numerical algorithms ifTTl . lfT2ll to compute Ci,j, here 
we derive closed-form j expressions for K up to four with 
arbitrary N . Note that the considered cases K < A cover the 
typical situations in applications discussed in Section |I] 

Following the methodology of obtaining the coefficients for 
the smallest eigenvalue distribution 1.14,1 . we first write an 
integral representation for the largest eigenvalue distribution; 

D{K,N) 



{K-l)\ 

K 



n i^^-^jf 



2<i<j<K 



N-K, 



{x - XifdX^, 



(27) 



i=2 



Cl. 



K = 2: 
2{~IY-^{2N - j - 2)\ 
{-N + j + 2)\{N - j)\{N - 2)\{N - 1)! ' 



C2. 



-{2N-]-2){2N-j-i) 
{-N + 3 + 2)\{N ■ 



(23) 



C2, 



C3, 



iv: = 3 : 

""^'v^^'''^ 2(-lp-^+HA^ - fc)!(2jV ^3 + k- A)\{-N + j - 2fc + 4) 

fe=max{^-jv+i} (2 - - J + fc - +3-k + 3)!(7V - 3)!(7V - 2)!(iV - 1)! ' 

(-1)^ / 2(jV-fc)!(2jV-j + fc~3)! 

^ (2 - A:)!fc!(-iV + ? - fc + 3)! I (2A^ - ? + /c - 5)!(A^ - 3)!(iV - 1)! 

(iV-fc- l)!(2iV-j + fc-2)! (iV-fc + l)!(2A^-j + /c-4)! \ 



(2iV - j + fc - 4)!(iV - 3)!(Ar - 2)! (27V - 3 + k- 6)!(7V - 2)!(Ar - 1)! / ' 

mm{j-N+3,N-l} 



E 



1 / (A^-fc + l)!(2iV- j + fc-4)! 



2fc!(-Ar + 7 - /c + 3)!(A^- 2)! \ (A^ - - 1)!(27V - i + fc - 6)! 
(N - ky.{2N - j + k - 3y.{N - 2) \ 



{N -k- 2)!(2iV -j + k- 5)!(iV - 1) 



(25) 



(26) 



where the constant D{K, N) = [l\f^i{N - i)\{K - i)\ ' 

and the domain of the integration J — [0,a;]^^^. Similar to 
the case of the smallest eigenvalue lfT4l . we first define the 
following integral 



La{x) := / X^ix-Xfc-^dX, 
Jo 



(28) 



which, by repeated use of integration by parts, equals 

' 2(a-fc + 2)! , 



L 



(-l)'^-fc!(2-fc)!' 



(a-fc + 2)!^,^ 



kl{a-ky. 



(29) 



fe=0 ^ ' ' k=Q 

When K ^ 2, the distribution in ( p7] i becomes 

fx, (x) = i?(2, iV)x^-2e-"iw_2(x). (30) 

Comparing ( [30| with (j6]l and after some manipulations, the 
coefficients are obtained as ( 23 1 on top of this page. 



For K = 3, it can be verified that ( 27 1 equals 

h,{x) ^ D{3,N)x^-''e'^(^LN-i{x)LN-3{x) 

iLN-2{x)f 

By using the equality 

a b a+b min{i,a} 

Ep«^'E^*^*"E E Pkqi~kX\ 

1=0 1=0 fc=max{0,i-6} 



(31) 



(32) 



i=0 



and comparing ( 3 1 1 with d6]i the coefficient expressions can 



be calculated as shown in p4]i-( 26 1 on top of this page 



For K — A, equation ( 27 1 can now be represented as 

fxAx) - I)(4,iV)x^-V"(^2Ljv-i(x)L^_2(:j)iA'-3(2:)- 
LN{x)LN-2{x)LN-i{x) - {Ln-3{x))'^Ln{x) - 



{Ln-i{x))^Ln-4:{x) - {LN-2{x)y 



(33) 



Using the equality 



'^p,x'^q^x'^kx' = 

1=0 1=0 1=0 

a+b+c mm{i,a+b} min{t,a} 

E E E 

2=0 t=max{0,i— c} fc=max{0,t — 6} 



Pkqt-kh-tx\ (34) 



the closed-form coefficients for K = A can be similarly 
obtained. They are, however, omitted in this paper due to space 
limitations. 

Note that for K >'b the numerical algorithm outlined fTT], 
l[T2| needs to be used to obtain the coefficients. Interested 
readers may contact the first author for a copy of the code of 
the numerical algorithm implemented in Mathematica®. 

IV. Numerical Results 

Some numerical examples are provided in this section to 
verify the derived SLE PDF expression ( [T3| ), CDF expres- 
sion ( [T4| and corresponding coefficient expressions. We first 
examine the cases when both closed-form distribution and 
coefficients are available, where we choose iV = 10 with 
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Fig. 2. CDF of the scaled largest eigenvalue using closed-form coefficients. 




Fig. 3. PDF and CDF of the scaled largest eigenvalue: K 
using numerical coefficients. 



N = &, 



K = 2,3,4 respectively. Using the derived CDF expres- 
sion ( [T4| with the corresponding closed-form coefficients, the 
CDFs against simulations are plotted in Fig. [2] In addition, 
we consider the case when K = Q and = 6, where 
the closed-form coefficients are not available. In this case, 
inserting these numerically obtained coefficients from ifTTIl 
(Table IV) into ( pj) and ( [T4| , the corresponding PDF and 
CDF are drawn in Fig. [3] From Fig. [2] and Fig. [3] we can see 
that the derived results match simulations well. 



results are easy and efficient to compute, they do not involve 
any complicated integral representations or unexplicit expres- 
sions as opposed to existing results. The derived expressions 
were validated through simulation results. 
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V. Conclusions 

Knowledge on the statistical property of the scaled largest 
eigenvalue of Wishart matrices is key to understanding the 
performance of various hypothesis testing procedures and 
communication systems. In this work, we derived exact ex- 
pressions for the PDF and CDF of the SLE for arbitrary matrix 
dimensions by using a Mellin transform based method. Our 



